Theory of the Magnetic Moment in Iron Pnictides 
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We show that the combined effects of spin-orbit, monoclinic distortion, and p-d hybridization in 
tetrahedrally coordinated Fe in LaFeAsO invalidates the naive Hund's rule filling of the Fe d-levels. 
The two highest occupied levels have one electron each but as a result of differing p-d hybridizations, 
the upper level is more itinerant while electrons in the lower level are more localized. The resulting 
magnetic moment is highly anisotropic with an in-plane value of 0.25 — 0.35ps per Fe and a z- 
projection of 0.06/is, both of which are in agreement with experiment. 

PACS numbers: 71.10Hf,71.55.-i,75.20.Hr,71.27.-|-a 



The representative parent material, LaFeAsO, in the 
rapidly growing class of iron-based superconductorsii^i^ 
exhibits a structural monoclinic distortion from tetrag- 
onal symmetry at IhQK followed by a transition to 
an antiferromagnelj^'i'^ at 134K with a unit cell of 
(%/2 X V2 X 2). The observed magnetic moment per 
Fe atom has been reported to range from 0.25/iB^ to 
and lies in the a-b plane. Such a low value of the 
magnetic moment is astounding because any application 
of Hund's rule to the Fe d-states results in a moment of at 
least 2hb ■ We offer here a resolution of the low in-plane 
magnetic moment in LaFeAsO which is rooted in three 
effects that are well known to be important in FeAs-based 
materials^iiSiiii, namely spin-orbit, strong hybridization 
between the Fe d and the As 4p orbitals, and the lat- 
tice compression along the z-axis to the lower monoclinic 
symmetry. All three conspire to destroy the naive Hund's 
rule filling of the Fe atomic levels as illustrated in Fig. [T] 

That the properties of Fe-based materials are strong 
functions of the hybridization is not new. Well known is 
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FIG. 1: Evolution of the energy levels of the Fe 3d and As 4p 
levels after the inclusion of spin-orbit coupling, p-d hybridis- 
ation, Vpd, and the monoclinic crystal field distortion. 



the case of an isolated Fe atom which has a moment of 
4/iB whereas in the metal, the moment is roughly halved 
to 2.2/iB per Fe as a resulti? of the 4s and 3d hybridiza- 
tion. Less well-known, but more pertinent to LaFeAsO, is 
the case of the Zincblende complex FeAs which is also an 
antiferromagnet and has a monoclinic distortioni^. First 
principles calculations on Fe films deposited on GaAsiS 
reveal that the value of the magnetic moment per Fe is 
a strong function of Fe-As bond distance. The moment 
vanishes^ for Fe-As distances less than 2.36A. This ef- 
fect was attributedi^ to the strong hybridization between 
the Fe 3d and the As 4p orbitals. In LaFeAsO, the av- 
erage Fe-As bond distance, 2.4A, is close to the critical 
value of 2.36A found for Fe-As films. As the degree to 
which Fe and As are non-coplanar in FeAs and LaoFeAs 
is identical, similar extreme sensitivity of the moment to 
the p-d hybridization is expected in LaFeAsO. 

It is ultimately symmetry that dictates hybridization. 
In LaFeAsO, each Fe is tetrahedrally coordinated. Full 
tetrahedral (cubic) symmetry splits the d-states into two 
irreducible representations^ili: 1) the three- fold degen- 
erate Fis levels consisting of the dxy, dyz, and dxz and 
2) the doubly degenerate consisting of dx2_y2 and 
dz2. The levels lie lower in energy. It is important to 
note that in a tetrahedral field, only the Fis states have 
the right symmetry to hybridize with the p-states of the 
sp neighbouring atom (As in this case), forming bond- 
ing and antibonding hybrid orbitals. The levels re- 
main non-bonding and hence will be neglected in our hy- 
bridization analysis. They will be assumed to constitute 
a full band (4 electrons). The immediate problem with 
applying Hund's rule to the Fis levels is that the effective 
moment on these levels is at least 2yUs as found in recent 
calculations^^. While inclusion of magnetic frustratiorJ^ 
might offer some reduction in the moment, it offers no 
resolution of the problem that the moment lies in the xy 
plane. The answer lies elsewhere as suggested by recent 
first- principles calculations^ - and a p-d mixing model^^. 

The experimental observation that the Fe moment lies 
in the plane is highly suggestive of spin-orbit coupling. 
To this end, our starting point is a general model, 

H = pV(2m) + Vo + n/{Arr?c^) {VVq x p) • S , (1) 
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for a cubic crystal with spin-orbit interaction where p is 
the momentum operator and S is the spin operator. This 
interaction breaks SU{2) symmetry. The rough idea is 
to include the effects of p-d hybridization and the z-axis 
distortion through a series of successive diagonalizations 
to obtain the eigenstates in the final basis. We only in- 
clude the outline of this calculation since an analogous 
analysis has been done for chalcopyrite semiconductors^. 
In obtaining the basis that diagonalizes the spin orbit in- 
teraction, we define |±) = \{X ±iY) / ^2) ,\0) = \Z) which 
are eigenstates of angular momentum {L, L^) with eigen- 
values (1,±1) and (1,0) respectively. The Hamiltonian 
for the p-levels is diagonalised through 
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States with the same indices are degenerate and r„ de- 
notes the symmetry of a state. Likewise, the basis for 
the d-levels which initially have Fis symmetry. 
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is formed from the states, |©) = \{YZ + iZX)/V2), 
|0) = \{YZ - iZX)/V2), |0) = \XY) are the eigen- 
states of {L,Lz) with eigenvalues of (2, +1),(2, — 1) and 
(2,0), respectively. As is clear, none of these states is 
an eigenstate of as is expected once the SU(2) spin 
symmetry is broken by the spin-orbit interaction. 

To consider the hybridization, we collate the states into 
two groups, segregating them according to their super- 
script a or f3. Within each group they are ordered as 
follows: 0^i(r8),02i(r8), 0^o(r8),0So(r8), ^^rr) and 
(f>22(^7)- Taking into consideration that only states of 
the same symmetry mix, 

(0p.(r™)|^pdlC,(r„)) = AM.,<5 (4) 

we find that the hybridization matrix can be written as. 
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Ap and Ag are the spin-orbit splitting of the p and d 
band, respectively. The highest three eigen-energies are, 

Ao = Ai = [Ap/3 + E- Arf/3] /2 + ^^7/2, (6) 



A2 = [-2Ap/3 + + 2A<j/3] /2 + V^/2 , 



(7) 



where ai = (A/3-i; + Ad/3)2-h4M2 and a2 = (2A/3 + 
-E+2Ad/3)^+4M^. According to the symmetries r„, n — 
7 or n = 8, the corresponding eigenstates are <I>f (r„) = 
ai0pj(r„) + 6i(/)^-(r„) with the coefficients {i = 0, 1, 2) 
and bi defined as. 



70 = 71 = a? = 1 - = [1 + MV {\^-E + Ad/3)' 
72 = 02 = 1 - ^2 = [1 + M'^l {^2-E- 2Ad/3)' 



To gain information about the spins, we transform the 
operator for the z-projection of the spin. 
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and cr/ = — cr^" into the or Ip basis. 

Ultimately, we will focus only on the three highest 
eigenstates. We refer to this reduced basis as II"'''. The 
final ingredient is the z-axis distortion from perfect cubic 
symmetry. We consider a crystal field interaction of the 
form, 

{X\V,i\X) = C{Y\Va\Y) = 5p/3, {Z\V,i\Z) = -2,5p/3 , 

{zx\v,i\zx) = (:{yz\v,i\YZ)^5d/i, 

{XY\V,i\XY) = -2(5p/3. 

The parameter C, accounts for the distortion in the a — h 
plane. Experimentally, the lattice constants along a and h 
differ by 0.3%. While any distortion is sufficient to lower 
the U{1) rotational symmetry in the plane to simply Z2 
(Ising), this effect is small relative to the overall z-axis 
tilt. As the parameter 5p is certainly not known within 
0.3%, we consider onlythe case of C = 1. The crystal field 
Hamiltonian in the 11^''^ is, 



T4f = 



Ai - Ti T2 
Ao+Ti 
T2 A2 



(9) 



whe re Fi = \ [5p-ii+ 5d{l - 71)] and T2 = ^[5py/rn2 + 
6d\/ (1 — 7i)(l — 72)]- Diagonalizing this Hamiltonian 
gives rise the following three energy levels, Eo{Tq) = 
Ao + Fi and 
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and their corresponding eigenstates, 
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which we refer to as the III 



basis, where i = 1,2 and 



Ci and are defined as, = 1-dl = l/[l+r|/(i?i-A2)] 
with (^2 = — ci, and C2 = di- In the final basis, III"*-'^-', 
the cTz spin matrix becomes 



TABLE I: Energy levels, hopping matrix elements, and in- 
teractions in the three highest levels in transformed basis, 
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where e = 2cidir + (df - c?)/3, / = (d? - cl)T - 2cidi/3 
and r = 2-\/2(aia2 + &i&2)/3. So we can see clearly that 
the final basis does not diagonalise the the spin matrix. 
Consequently, the states 111"^'^^ represent some linear 
combination of spin up and spin down. A crucial point 
about this spin matrix: the a and (i states have opposite 
projections of spin in the states with energies Ei and E2- 
In the transformed basis, the resultant Hamiltonian 
reads. 
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uncertainty in the value of these parameters, especially 
the hybridization energy, M and the monoclinic distor- 
tion. Hence, we explore the dependence of the Hamilto- 
nian parameters on both. As shown in the first panel of 
Fig. IHa), the E2 level is the lowest followed by Eq and 
then El for sufficiently large values of M . As can be seen 
from Fig. mb), U2 > E0 — E2 and Ei—E2- Consequently, 
E2 and one of Eq or Ei will be singly occupied. To deter- 
mine the ground state spin configuration, we note that 



both the interactions U- 



a0 
12 



^ and Uo2^ 



^ (see Fig. 



(121;)) with a ^ for the highest two occupied levels. Re- 
call, the z-projection of the spins in a and /3 is reversed in 
levels Eq (or E\ ) relative to E2 ■ Level Eq is an eigenstate 
of Sz while El is not. Consequently, the lowest energy 
configuration for single occupancy of the levels Eq and E2 
is an antiparallel alignment of the spins. That is, both 
electrons are in the a or in the (3 state in both levels. 
This configuration does not cost the repulsion term Uq2 
or U12 ■ Since there are 12 electrons to fill these levels, 
we arrive at the filling structure shown in Fig. [TJ Con- 
sequently, the three effects considered preclude a naive 
assignment of the spins according to Hund's rule--^. This 
conclusion is a general result of this analysis, not an ar- 
tifact of fine-tuning the bare parameters. 



All other couplings, for example, toi and ^02 vanish by 
symmetry. The coefficients Ai,Bi,Ci and Di are defined 
as. 
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\/2j3c.iai 



\f2fidia1, 
\l\fZdiai, 



Di 



1/V3c,bi - 
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where ci, di, a;, hi are defined as before. 

We analyzed the energy levels, interaction strengths, 
and spin projections based on representative values for 
iron-based systems. For instance, if we set Up = Ud = 
4eF.i^, tp = td = O.TeV-^^, M = O.SeV^^, Ap = 
0A26eV^, Ad = O.OSeV-i^, E ^ 1.9eV^, and 6p = 
5d = O.OGeV"— we arrive at the parameters of Table I 
for Hamiltonian Notice, however, that there is an 



The problem has now been reduced to the physics of 
two low- lying energy levels, Ei and E2 ■ That a two-band 
reduction reproduces^i the Fermi surface seen in the local 
: 2) density approximation^^ and experiments^'^ corroborates 

2/1 our approach. If we neglect the interactions in (jl3p . the 
\B1B2) /"AUt^qY^Xquy can be easily diagonalized and one finds two, 
(ld£)ubly degenerate, energy-shifted bands. The sphtting 
between the bands is due to the crystal field that shifts 
each half-filled band away from the perfect nesting con- 
dition. As a result, it is difficult to reconcile the exper- 
imental observation of a spin-density wave (SDW) with 

a simple Fermi surface instability due to nesting. The 

V 2/3<ii&i, resolution may lie in the interactions. As Fig. [D(b) in- 
l/\/idibi, dicates, for all values of the hybridization, the levels 1 
and 2 have differing p and d character. In E2 the on- 
site interaction, U2 exceeds Uq (or Ui) by more than a 
factor of two at M = O.SeV. Consequently, single oc- 
cupancy in the Eq (or the Ei) level results in itineracy 
whereas in the E2 level Mott physics can be relevant since 
U2 ~ 4^2. This difference is due entirely to the different 
p-d character between the Eq (or Ei ) and E2 levels which 
is expected as they are orthogonal. Finally, we show in 
the last panel in Fig. [HJd) the value of Sz in the lev- 
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FIG. 2: Energy levels, hoping matrix elements, on-site in- 
teractions, intra-band interaction and spin components as a 
function of hybridization M. Shown in the inset of (c) is the 
on-site interaction U2 as a function of the monochnic distor- 
tion. 



els El and £2- Recall level Eq is an eigenstate of Sz 
with the z— projection opposite to that in i?2. As shown 
in the last panel in Fig. ([2]), the value of Sz in E2 is 
If Eq is the next lowest level, a net moment in 
the z— direction of remains as has been recently 

observed^. If Ex is relevant (as would be the case for 
M > 0.83eF), the z-moment vanishes as shown in Fig. 

. The itineracy of the electrons in the Eq level does 
not affect this cancellation as it is the average not the lo- 
cal spin configuration that is relevant in a magnetization 



measurement. Hence, for experimentally relevant values 
of M {.heV < M < l.OeV), the residual z-component 
of the moment is strongly dimini shed. Th e remaining 
x-y component on E2 is Sxy = y'l — {Sz}"^ — 0.317/is. 
Such a moment can order via the super-exchange mech- 
anism on E2 as U2 ^ ^2- The evolution of S^y as a 
function of A/ in Fig.[2ljd) shows that the analysis here is 
consistent with the range of the magnetic moment seen 
experimentally^. 

Our analysis also explains why the structural 
transition^ at 150K must precede the onset of antifer- 
romagnetism. As the inset in Fig. [2Uc) indicates, the on- 
site energy U2 diminishes as the crystal field associated 
with the monoclinic distortion vanishes. Once U2 < t2, sl 
transition to an antiferromagnet via the super-exchange 
interaction is untenable. The structural transition breaks 
rotational symmetry in the plane not SU(2) which is al- 
ready broken at the outset by spin orbit interaction. The 
success of the analysis presented here in describing the 
antiferromagnet in the parent material, LaFeAsO, im- 
plies that Eq. p3p should be used in any subsequent 
analysis of superconductivity. The presence of an itin- 
erant band coupled to one with moderate Mott physics 
makes the problem of the iron pnictides more akir^ to 
that of the Kondo lattice in heavy fermions than the 
cuprates. 
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Note Added While this work was under review, 
McGuire, et al^ reported a residual magnetic moment 
along the c-axis equal to 0.06/z_b as predicted here. 
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